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Abstract This experimental study reports on the be¬ 
havior of a filament attached to the rear of a three- 
dimensional cylinder. The axis of the cylinder is placed 
normal to a uniform incoming flow and the filament is 
free to move in the cylinder wake. The mean position 
of the filament is studied as a function of the filament 
length L. It is found that for long ( L/D > 6.5, where D 
is the cylinder diameter) and short ( L/D < 2) filaments 
the mean position of the filament tends to align with 
the incoming flow, whereas for intermediate filament 
lengths (2 < L/D < 6.5) the filament lies down on the 
cylinder and tends to align with the cylinder axis. The 
underlying mechanism of the bifurcations are discussed 
and related to buckling and inverted-pendulum-like in¬ 
stabilities. 

Keywords Wake • Symmetry breaking • Fluid- 
Structure interaction 


1 Introduction 

Many animals, plants and seeds have appendages or 
protrusions that through a passive interaction with the 
surrounding fluid can aid locomotion. The underlying 
physical mechanisms responsible for such positive fluid- 
structure interaction can be isolated using simple body- 


flow, it is found that for long lengths of the splitter 
plate, the mean position of the splitter plate is sym¬ 
metric - i.e. it aligns with the incoming flow - whereas 
for short lengths, the symmetry is broken, i.e. the mean 
position of the splitter plate will form an ang l e with 
the incoming flow. More recently, iBagheri et all ( 20121) 
studied through two-dimensional numerical simulations 
the behavior of a single elastic filament hinged on the 
rear of a fixed bluff body and placed in a uniform in¬ 
coming flow. They also found a symmetry-breaking for 
sufficiently short filament lengths and reported on the 
induced lift/side force acting on the cylinder-filament 
system. More specifically, they showed that for a long 
filament (symmetric position), the mean lift force and 
torque acting on the cylinder is zero, whereas when the 
filament length is reduced below two cylinder diame¬ 
ters, a significant mean lift force/torque ac ting on the 
cylind er is generated. Even more recently, Lacis et a3 
( 20141) performed numerical simulations of a free-falling 
circular cylinder with a splitter plate and found that 
for short splitter plates the body drifts - as a conse¬ 
quence of symmetry breaking - either to the left or 
right. They also proposed a theoretical model - referred 
to as the inverted-pendulum-like (IPL) model - that 
provides quantitative predictions of the critical filament 
length of the symmetry breaking. In this model, the 
pressure in the recirculation zone behind the bluff body 
acts n s a destabilizi ng force and is res ponsible for the 


appendage configurations . For example, two- dimensional 

nume rical (|Xu et al|,[l99o|) and experimental dCimbala and Gjheiahtry breaking. lLacis et all d2014h also found nu- 
19941) studies have reported on a cylinder free to rotate 
around its axis with a rigid splitter plate attached to 
it. When this system is placed in a uniform incoming 
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merically the IPL instability in a three-dimensional con¬ 
figuration consisting of a sphere (free to rotate and 
drift) with an attached thin elliptic-shaped sheet. It 
was confirmed that the mean position of the short sheet 
forms an angle with the incoming flow and thus drifts. 
In summary, previous studies have shown that flexi¬ 
ble/rigid appendages can be used to passively control 
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the flow, for example by mod ifying the trajectory of 
fallin g bodies such as cyli nders (iHorowitz and Williamsonl . 
201 oh , spheres or disks ( Ern et al . 2012h . The aim of the 
present study is to experimentally investigate the be¬ 
havior of a three-dimensional cylinder with a filament 
attached to its rear and placed in a uniform flow. To the 
best knowledge of the authors, this configuration has 
not been investigated previously, as earlier work has 
either been two-dimensional or numerical. The physi¬ 
cal problem depends on seven dimensional parameters: 
the cylinder diameter D , the length of the filament L, 
the fluid density pf and filament density p s , the bend¬ 
ing stiffness of the filament B , the free-stream velocity 
U and the kinematic viscosity v. Four non-dimensional 
parameters may then be introduced, namely, the length 
ratio L s = L/D, the density ratio p = p s /pf , the 
Reynolds number Re = UD/u and the bending ra¬ 
tio R\ = B / (pfU 2 D 4 ). In the present study, two free- 
stream velocities (U = 0.92 cm.s -1 and U = 1.56 cm.s -1 ) 
are investigated, which results in two configurations 
(Re = 185, R 1 « 7.4 x 10 -4 ) and (Re = 310, R 1 « 
2.6 x 10 -4 ). For each of these configurations, the length 
ratio is varied 0.5 < L s < 15, but the density ratio is 
kept fixed (p = 1.25). 


2 Method 


The experiments have been carried out in a vertical wa¬ 
ter tunnel; a schematic of the experiment is presented 
in figure [U The whole apparatus is about 1.5 m high, 
1 m wide and 0.4 m deep. The dimensions of the test 
section are 20.5 cm wide, 39.7 cm deep, 68 cm high. Be¬ 
low the test section a 40 cm long slab of foam is placed 
in order to homogenize the flow. The vertical flow of 
water is generated with a pump with typical flow rate 
of0.76L.s _1 , which corresponds to an upward fluid ve¬ 
locity U = 1 cm.s -1 in the test section. The cylinder is 
made from a plexiglass tube with a diameter D — 20 
mm and length of 20.5 cm. This cylinder is placed in 
the test section with its axis of symmetry normal to the 
incoming flow. A hole is drilled in the cylinder through 
which a thread made from silk is inserted. The flow 
inside the channel was previously measured by iBellanil 
( 2005 1 with ultrasonic velocity profiler (UVP). At the 
position of the cylinder - 20 cm above the foam - fluid 
velocity profiles for U z were measured along the y di¬ 
rection; the fluid velocity U z was found to vary within 
±5% over 50% of the center of the channel, the turbu¬ 
lence intensity of U z was found to be 5%. 

The position of the filament can be defined by two 
angles of deflection, one from the front and one from 
the side - called 0\ and # 2 , respectively. These two an¬ 
gles are shown schematically in figures Eta) and Eta). 



Fig. 1 Schematic of the vertical water tunnel facility. 


As the shape of the filament is not a straight line but 
can be curved it can be noted that the two angles of 
deflection 6 \ and 62 do not take into account the shape 
of the filament. As a consequence, the two angles give 
only partial information regarding the position of the 
filament, nevertheless 6 \ and 62 indicate the behav¬ 
ior of the filament regarding the symmetry breaking 
and the lie down behavior that this paper focuses on. 
In order to measure the position of the filament, the 
cylinder is photographed with two perpendicular cam¬ 
eras (1400 x 1088 pixels), one from the front and one 
from the side. The conversion factor is found to be 5 
pixels per mm corresponding to an accuracy of ±0.1 
mm in the determination of the filament position. For 
each length of the filament a series of grey scale im¬ 
ages are recorded; in total seventeen filament lengths 
are recorded between 0.5 < L s < 15. For each series 
an acquiring rate of 2 Hz and an exposure time of 10 
ms was used; for the Reynolds number Re = 185 the 
recording time is 4 minutes and 10 seconds resulting 
in 500 images and for Re = 310 the recording time is 
8 minutes and 20 seconds resulting in 1000 images for 
each measured filament length. The greyscale images 
are then converted to black and white images in order 
to detect the filament position. The angles of deflection 
are calculated as the angles between the vertical and 
the line formed by the filament attachment point to 
the cylinder and the endpoint of the filament (positive 
angles are defined in the clockwise direction). To obtain 
the mean angles of deflection, the angles for the entire 
measurement series (500 and 1000 images for Re = 185 
and Re = 310 respectively) are averaged in time and 


























Experimental study of a three dimensional cylinder-filament system 


3 



Fig. 2 (a) Schematic of the cylinder-filament for the front 

view The grey region marks the back-flow region. Inside this 
region the flow is in reverse direction to the free stream. The 
normal forces on the filament inside and outside the 
back-flow region and the tangential forces F+ and F t - are 
also shown; (b) photography of cylinder-filament system for 
the front view. 




Fig. 3 (a) Schematic of the cylinder-filament for the side 

view; (b) photography of cylinder-filament system for the side 
view. 

the absolute value of this mean is then taken to give 
the mean angles of deflection (9 1 ) and (# 2 )- 

3 Filament behavior 

In this section, the time evolution of the filament po¬ 
sition in the cylinder wake are shown, followed by a 
study of the mean position of the filament as a func¬ 
tion of its length. Figures 0J [5] and [ 6 ] show the time 
evolution of the angles of deflection, 9\ and 62 and the 
corresponding filament positions at different times for 
three filament lengths, namely, L s = 9.2, L s = 4.1 and 
L s = 2.1 (Re = 310). For the longest length L s = 9.2 
shown in figure 0J the angles 9\ and O 2 oscillate around 
a value close to 0 °, indicating that the filament is al¬ 


most vertical (figures 0^ and 03 ). This can also be seen 
in the filament positions presented in figure [U^c) for the 
front view and in figure 0Jd) for the side view. Fig¬ 
ure [5] shows that for the intermediate filament length, 
L s =4.1, the angle 9\ oscillates around a value close to 
23° and the angle 62 oscillates around 80° (figures [5^ 
and [Sb). This indicates that the filament is lying down 
as can be observed in figures [5jc) and[5fd). Note that 
for this situation, the detectable part of the filament 
from the front view is reduced; however, the angle 9 2 
clearly shows the collapse of the filament towards the 
cylinder surface. For the shortest filament L s = 2.1, the 
angles 9\ and #2 are close to 0 ° (figures [ 6 ^ and [ 6 ) 3 ) in¬ 
dicating that the filament is again almost vertical. This 
can also be observed in the filament position in figure 
[ 6 ](c) and[ 6 ](d). 

The average angles of deflection (9 1 ) and (# 2 ) are 
shown in figures 0a) and 0b) as a function of L s for 
the two Reynolds numbers under investigation. For fil¬ 
aments longer than L s = 6.5, (# 1 ) and ( 62 ) are smaller 
than 5° and 25°, respectively. This indicates that the 
mean position of the filament is nearly vertical from the 
front view - i.e. aligned with the incoming flow - and 
slightly tilted from the side view (schematically shown 
as inset frames in figure 0 . Figure [ 8 ] presents the stan¬ 
dard deviation of the angles 9\ and #2 as a function of 
the filament length for Re = 310. For these long fil¬ 
ament lengths the standard deviation of the filament 
oscillations are below 5° for both angles and increase 
close to the bifurcation. 

When the filament length is reduced below a criti¬ 
cal value, L s ~ 6.5, a first bifurcation takes place as the 
angle (# 2 ) increases dramatically to reach values higher 
than 70°. The angle (# 1 ) also increases at this bifurca¬ 
tion. Thus, in this regime the filament tends to collapse 
(i.e. fall down) on the cylinder and align with the axis 
of symmetry of the cylinder (normal to the incoming 
flow) as schematically shown in figure 0 Once the fila¬ 
ment collapses on one side, it does not flip to the other 
side. Out of 10 experiments we performed for this par¬ 
ticular configuration, the filament fell down to the left 
twice and to the right 8 times. In perfect symmetric 
conditions, one would expect an equal probability for 
the filament to fall to right and left, however further 
samples and statistical analysis are required to confirm 
this. For these intermediary filament lengths the stan¬ 
dard deviation of the angle 9\ varies between 5° and 
21 °; the standard deviation of the angle #2 varies be¬ 
tween 0.5° and 13° (figure [ 8 ]). The evolution of these 
amplitudes with the filament length are non-monotonic 
in this range. 

Finally a second bifurcation takes place at L s ~ 3 
for Re = 310 and at L s ~ 1.5 for Re = 185; the value of 
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Fig. 4 The angles of deflection for (a) 0 1 and (b) 62 as a 
function of time and filament position at different times {t = 0 
s to t = 10 s) for (c) the front view and (d) the side view; 
L/D = 9.2 and Re = 310. 



Fig. 5 Same as figure 2| but for L/D = 4.1. 


the mean angle (0 1 ) is smaller than 20° (figure 0i) and 
the angle ( 62 ) is reduced from values higher than 65° 
to values lower than 20° (figure 03). This indicates that 
the mean position of the filament is still tilted, but tends 
to realign with the incoming flow as shown schemati¬ 
cally in figure 0 For these short filament lengths the 
standard deviation of 6fi and O 2 decrease from respec¬ 
tively 7° and 25° to 0.5° for both angles (figure [8]), in¬ 
dicating that the amplitudes of the filament oscillations 
decrease with the filament length. 

Concerning the frequencies of the filament oscilla¬ 
tions, it can be observed from figures 0(a), 0(b), 0(a), 
0(b),[61(a) and 0(b) that several distinct frequencies are 
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x/D y/D 

Fig. 6 Same as figure [4] but for L/D = 2.1. 


contained in each signal. For 62 no clear frequency seem 
to dominate. For on the other hand, the Strouhal 
number St = fiD/U of the dominant frequency /1 is 
plotted as a function of the filament length in figure [9] 
(for Re = 310). For long filament lengths, L/D > 6.5, 
the Strouhal number of is between 0.14 and 0.2, 
which is close to the natural vortex shedding frequency 
of the cylinder wak e St cl , under ~ 0.19 (at Re = 300, 
see e.g. fwilliamsonl (|l996h ). Thus, for L/D > 6.5 the 
filament oscillation frequency seems to be somewhat 


synchronized with the vortex shedding in the cylinder 
wake. Note however that contains also other frequen¬ 
cies that are not shown here. For intermediary filament 
lengths 2 < L/D < 6.5 no clear dominant frequency 
could be identified. Finally, for short filament length 
L/D < 2 the Strouhal number of the filament oscil¬ 
lation frequency is between 0.1 and 0.15 and is signifi¬ 
cantly different from St cy u n der • The different dynamical 
regimes determined from the average angles seem thus 
to be reflected in the dominant time scale of 6fl. 

We can also note some differences in the bifurca¬ 
tions between the two Reynolds numbers studied here, 
namely Re = 185 and Re = 310. As previously indi¬ 
cated, the main marker of the bifurcations is the angle 
62 and it can be compared between the two Reynolds 
numbers on figure 0(b). For the higher Reynolds num¬ 
ber, Re = 310, the first bifurcation appears at L s ~ 6.5 
and at L s ~ 3 for the second whereas for Re = 185 
the bifurcations appear at a slightly different filament 
lengths; at L s ~ 6 and L s ~ 1.5. The bifurcations 
seem also to be more sharp at Re = 310 compared 
to Re — 185 where it is more continuous. Despite these 
differences, it is clear that the same bifurcations are 
observed for both Reynolds numbers and these bifur¬ 
cations seem to be robust in this range of Reynolds 
number. 


4 Discussion 


Although, further quantitative investigations are nec¬ 
essary for confirmation, we believe that it is likely that 
the bifurcations described above are related to a buck¬ 
ling instability and an IPL instability. At the Reynolds 
numb ers considered in t his study, the recirculation re¬ 
gion ([Williamson . 1996 ) can be characterized by its 
mean length, which decreases from 1.5 to 1 D when the 
Reynolds number increases from Re — 120 to Re = 300 
(iSilva et all . 2003 : Nishioka and Sato . 1978 : Park et al 
19981 ). Inside this recirculation region we can define a 
back-flow region where the fluid velocity is opposite 
to the streamwise direction. This back-flo w region can 
be m odelled as an ellipse in the xz-plane ( Lacis et~al 
20141 ) and in the three-dimensional case this ellipse is 
extended in the z—direction. The shape of the back-flow 
region is schematically shown in figure [2fa) for the xz- 
plane and in figure [31(a) for the yz- plane. The behaviour 
of the filament can be related to this mean flow and may 
be attributed to two different types of instabilities. 

The first instability is a simi lar mechanism a s the 
two dimensional IPL instability ( Lacis et"al 2014 ). i.e. 
the pressure force in the back-flow region behind the 
cylinder, exerts a net normal force F+ on the filament, 
which pushes the filament away from the symmetric 
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vertical position. The part of the filament outside the 
back-flow-region on the other hand, experiences a net 
normal restoring force F~, which pushes the filament 
towards the straight vertical position. As shown schemat¬ 
ically in figure Ufa), the final position of the filament is 
due to the competition between these two forces, which 
for sufficiently short filaments results in a tilted mean 
position with respect to the directio n of the incoming 
flow. The key difference to the work of Lacis et a] ( 2014 ) 
is that in the present three-dimensional case, we have 
observed that the filament falls down on the cylinder 
in the vertical yz- plane containing the cylinder axis (in 
contrast to the xz-plane). In the yz- plane, the back-flow 
region has on average a rectangular form (in contrast to 
the elliptic shape in the xz- plane). As a consequence, 
the restoring normal force F~ outside the back-flow 
region becomes increasingly small as 62 increases and 
once the straight position of the filament is unstable 
no other force can prevent the filament to he down on 
the cylinder. Although, physically an IPL-type of in¬ 
stability could be responsible for the first bifurcation, 
the critical value of filament length near L s = 6.5 is 
much larger than the critical values (L s ~ 2) reported 
in previous work on the IPL instability and also much 
larger than the length of the recircula tion region which 
is ^ 1 D at these Reyno l ds numbers ( Silva et all . 2003 : 


Nishioka and Sat et Il978l : iPark et all . 11998 ) 


The second instability that could induce the col¬ 
lapse of the filament onto the cylinder surface in the 
xz-plane is a buckling-type of mechanism. This instabil¬ 
ity arises out of the competition between the tangential 
forces F t + and Ff (shown schematically in figure [3^) 
on the filament inside and outside the back-flow region. 
As a consequence of these forces, part of the filament 
outside the back-flow region is stretched by the flow 
and the part of the filament inside the back-flow region 
is compressed. If the filament is sufficiently short, the 
compressive force will dominate, which in turn may re¬ 
sult in a buckling instability of the filament. For the 
second bifurcation at L s ~ 2, where the filament tends 
to realign with the flow, the filament is unlikely to un¬ 
dergo a buckling instability as the critical load for the 
onset of the buckling instability scales as 1/L 2 . 


The above discussion suggests that the first bifurca¬ 
tion at L s = 6.5 could be dominated by a buckling-type 
of instability. Further studies, in particular by estimat¬ 
ing the loads acting on the filament, are necessary to 
confirm this statement and to quantitatively assess the 
mechanisms underlying the bifurcations. 




Fig. 7 Average angles of deflection (a) (0±) and (b) ( 62 ) as 
a function of the filament length for two Reynolds numbers 
Re = 185 and Re = 310. 
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Fig. 8 Standard deviation of the angles 6 \ and 62 as a func¬ 
tion of the filament length for the Re = 310 case. 
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Fig. 9 Strouhal number St = fijj of the oscillation fre¬ 
quency of Oi as a function of the filament length for the 
Re = 310 case compared to the Strouhal numbe r of the vor¬ 
tex s hedding in the cylinder wake at Re = 300 (IWilliamsonl . 
! 1996). 


5 Conclusions 

This experimental study has focused on the behavior of 
a filament attached to the rear of a three dimensional 
circular cylinder placed in a uniform incoming flow. The 
length of the filament was varied and for each length 
the mean position of the filament was determined. For 
lengths, L s longer than 6.5 the mean position of the 
filament is nearly aligned with the incoming flow. For 
intermediate lengths, 2 < L s < 6.5, it was found that 
the filament lies down on the cylinder, i.e. the mean 
position of the filament is aligned with the axis of the 
cylinder. Finally for L/D shorter than 2 the mean po¬ 
sition is tilted in both directions with respect to the 
incoming flow. To the best of our knowledge these two 
bifurcations of a flexible filament in wake flows have not 
been observed nor reported in earlier work. 
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